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ABSTRACT 

Photon-dominated regions (PDRs) are powerful molecular line emitters in external galaxies. 
They are expected in galaxies with high rates of massive star formation due to either starburst 
(SB) events or starburst coupled with active galactic nuclei (AGN) events. We have explored 
the PDR chemistry for a range of physical conditions representing a variety of galaxy types. 
Our main result is a demonstration of the sensitivity of the chemistry to changes in the physical 
conditions. We adopt crude estimates of relevant physical parameters for several galaxy types and 
use our models to predict suitable molecular tracers of those conditions. The set of recommended 
molecular tracers differs from that which we recommended for use in galaxies with embedded 
massive stars. Thus, molecular observations can in principle be used to distinguish between 
excitation by starburst and by SB+AGN in distant galaxies. Our recommendations are intended 
to be useful in preparing Herschel and ALMA proposals to identify sources of excitation in 
galaxies. 
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Galaxies may be powerful molecular line emit- 
ters because of a very high massive star formation 
rate induced either by starburst (SB) events or the 
presence of active galactic nucle i ( AGN) (jRadford 
Solomon et al.||1992| [Gerin fc Phillips 



subject of a modelling study ( Bayet et al.|2008b l. 
That work predicts that such galaxies should be 
strong emitters in rotational lines from molecules 
embedded in the warm cores from which the young 
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massive stars are evolving (see also Lintott et al. 
2005| |Lintott fc Viti|[2006| . |Bayet et al .| fl2008bl 



Bayet et al.|20 06 ; Evan s et al.|2006||Gracia-Carpio 
et al.||2006| I Martin et al.||2006[ |Aaito et al.||2007[ 



Bell et al. | |20"07l |Imanishi et al.| [20071 |Knudsen| 
et al.|2007|fPapadopoulos|2007||Krips et al.|2008| ). 



Some external sources may be dominated by only 
one type of activity or contain strong signatures of 
both. Both activity types (SB and AGN) are rec- 
ognized as having a direct influence on the molec- 
ular gas emission in a galaxy. 

Molecular emission from dense star-forming 
cores in starburst galaxies has recently been the 



described how the chemistry in these cores may 
depend on the local physical conditions, and they 
identified some molecular species that should be 
sensitive tracers of these conditions. 

Such starburst galaxies may also contain gi- 
ant photon-dominated regions (PDRs) which are 
also powerful emitters in molecular lines in the 
millimetre and submillimetre regions of the spec- 
trum. These PDRs are quite distinct physical and 
chemical systems from those discussed by [Bayet 
et al.| p008bl. Giant PDRs should be set up in 
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galaxies with clusters of young massive stars and 
in galaxies dominated by radiation from an ac- 
tive nucleus. It is therefore important to discuss 
the chemistry of galaxies dominated by PDRs, set 
up either by clusters of young massive stars or 
by the additional presence of active galactic nu- 
clei (AGN), and to identify the likely molecular 



tracers of such PDRs. In particular, it would be 
important to know whether these molecular signa- 
tures are similar to or distinct from those of the 
star-forming cores of starburst galaxies described 



by Bayet et al. (2008b). 



In this companion paper to Bayet et al. ( 2008b ) , 



we model the chemistry of PDR-dominated galax- 
ies. As described in the earlier paper, the physi- 
cal conditions that determine the chemistry in the 
interstellar medium of external galaxies (metallic- 
ity and/or relative elemental abundances; dust:gas 
ratio; dust grain properties; gas density; cosmic 
ray ionisation rate; and far UV radiation inten- 
sity) may differ substantially from those of the 
Milky Way and for our purposes are regarded as 
unknown parameters. The chemistry is explored 
for a wide range of these parameters and its sensi- 
tivity to variations in these parameters is studied. 
We also attempt to assign likely values of these 
parameters to galaxy type (starburst, SB+AGN 
and high-redshift) and predict suitable molecular 
tracers for these galaxy types. 

The PDR code and the chemical database are 
described in Section [2] and the physical param- 
eters and their ranges are discussed in Section 
[3] Results describing the PDR chemistry and its 
sensitivity to variations in the parameters listed 
above are presented in Section [4] In Section [5] we 
make predictions of the detectability of lines of cer- 
tain potential tracer molecules in PDR-dominated 
galaxies at low and high redshift, and in Section 
[6] we give a discussion of our results in relation to 
observations, to other PDR modelling, and to our 
earlier warm-core studies, and we summarize our 
conclusions. 

2. Model Description 

PDR models differ from each other in the as- 
sumptions adopted, in the geometry, or in the de- 
gree of sophistication used for solving the radia- 
tive transfer equations, or in the thermal balance 
and for determining the level population. PDR 
models also differ from each other by the chemi- 
cal network adopted and by the consideration (or 
not) of time-dependence. In this paper, we used 
the time-dependent UCL_PDR code benchmarked 
in Rollig et al. (20071 and fully described in Bell 



et al. (2006 20071 



dimensional semi-infinite slab is illuminated from 
one side by FUV photons. The chemistry and 
thermal balance arc calculated sclf-consistcntly 
at each point into the slab and at each time- 
step, producing chemical abundances, emission 
line strengths and gas temperatures as a func- 
tion of depth and time. The UCLJPDR code 
adopted here uses a chemical network contain- 
ing 131 species and over 1700 reactions from the 
UMIST99 database ( |Le Teuff et alpM)] ), includ- 
ing ion- molecule, photoionisation and photodisso- 
ciation reactions. The UCLJPDR model also in- 
cludes some modifications of UMIST99 introduced 
as part of the PDR benchmarking effort presented 
in Rollig et al. (20071. The freeze-out of atoms 
and molecules on to grains is neglected in this 
model. Apart from hydrogen (which is assumed 
to be initially mainly molecular), the gas is as- 
sumed to be initially in a purely atomic form, with 
initial elemental abundance ratios as free param- 
eters. Initial elemental abundances of all metals 
are assumed in most cases to scale linearly with 
the metallicity (z/z ), and the H2 formation rate 
and the adopted gas-to-dust mass ratio are also 
assumed to be linear with the metallicity. For a 
more detailed description of the UCLTTJR model, 



we refer the reader to Bell et al. ( 2006 2007 1 and 
references therein. 

3. Parameter choices 

To perform our study, we used and extended a 



model grid previously built and presented by Bell 
et al. (2006 20071. This grid contains more than 



1200 models, each with a different set of input pa- 
rameters. We restricted our selection to 9 models, 
the most relevant for our study. Table [1] lists the 
standard (Milky Way) values of the parameters. 
Table [2] lists the range of parameter values covered 
by the models we have selected. We adopt solar 
initial chemical abundances as our standard. Mod- 
els 9, 10 and 11 explore the effect of replacing the 
solar initial elemental abundance ratios by values 
predicted by early universe models from|Chieffi & 



The UCL_PDR model assumes that a one- 



Limongi|(|2002[);|Heger fc Woosley| ( |2002[ ); |Umeda 
fc Nomoto| Q2002| ) (see Table |3|). In these three 
models, the metallicity, the H 2 formation rate and 
the gas-to-dust mass ratio are fixed to standard 
values. As described later in Sect. [5] Model 12 is 
intended to represent high-redshift environments. 
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Models to 4 investigate the consequences of 
a reduction of solar metallicity by a factor up to 
100 consistently linked with a reduction of the H 2 
formation rate, of the gas-to-dust ratio and of the 
initial elemental abundance ratios. Models and 6 
study the consequent changes in the chemistry af- 
ter an increase of cosmic ray ionization rate from 
5 x lO-^s- 1 to 5 x 10- 15 s _1 . We suggest that 
the increase of cosmic ray ionization rate may be 
qualitatively used to simulate XDR-like environ- 
ments. Finally, a two orders of magnitude change 
in both the FUV radiation field (Models and 7) 
and in the gas density (Models and 8) have been 
investigated. 

Each model provides atomic and molecular 
abundances as a function of total cloud depth 
(at visual extinctions of < A„ < 10) for cloud 
ages of 10 4 , 10 5 , 10 6 , 10 7 and 10 8 yr. For the study 
of the A„ influence (Models and 5 in Table |2| , 
we have restricted our model selection to those 
showing either a value of A„=8, or a value of 
A^ = 3 (typical of a translucent gas component). 
Note that Model is considered as our reference 
model in the paper, assuming standard parame- 
ters (listed in Table [l| with an A„ = 8, likely to be 
a typical value for representing the dense PDR gas 
component in galaxies, as detected in the nuclei 



of M 82, NGC 253 , IC 342, and NGC 4038 ( |B-ayet 
e!TaT]|2008a[|2009l ). 



4. Sensitivity of chemical abundances to 
variations of the UCL_PDR model pa- 
rameters 

In this section, we present results on the trends 
in the chemistry with respect to the changes of var- 
ious parameters. We have selected 19 molecules 
for closer study, either for their likely detectabil- 
ity or for their chemical interest (or both). This 
set of 19 molecules is closely similar to the set of 
molecules used in Bayet et al. ( 2008b| ). We use 
this set to make a comparison between the predic- 
tions of the PDR model and those of models de- 
scribing the chemistry of high-mass star-forming 
regions (see Sect. |6|. We have arbitrarily fixed 
the limit of detectability of molecule X to be a rel- 
ative abundance of n(X)/n# = 1 x 10~ 12 (where 
n# =n(H) + 2n(Pi2)), a criterion roughly satis- 
fied in our own galaxy. The article deals specifi- 
cally with molecules CO, H 2 0, CN, OH, CS, HNC, 



HCN, HCO+, H 3 0+, SO, C 2 and C 2 H (see Figs. 
[I] to U and with C0 2 , OCS, S0 2 , H 2 S, H 2 CS, 
CH 2 CO and H 2 CO (see Tables [|-[§). Figs. [l]to[| 
show the time-evolution of molecular abundances 
and their sensitivity to parameter changes. Un- 
like for lower density models ( |Bell et al.|2006[ ), the 
chemical steady-state is achieved here in « 10 5 yr. 
Thus, time-dependence is unlikely to play a role 
in the chemistry. 

4.1. Metallicity 

Fig. [T] shows that the model with the lowest 
metallicity (Model 4 in Table [2| predicts very low 
abundances for all species except CO, H 2 0, HCO + 
and OH, and enhances significantly the abundance 
of H30 + . This enhancement is caused by a severe 
reduction in the electron density at low metallicity. 
Since the main loss of polyatomic ions is through 
dissociative recombination, the ion abundance in- 
creases. 

With z/z Q =0.01, we expect that the deuter- 
ated ions such as H 2 D + and H 2 DO + may be more 
readily detectable than H 3 + (see e.g. |Boreiko 



|fc Betz|[i"993] |Pagani et aT][2007| and references 
therein). The deuteration in both PDR and high- 
mass star forming models will be investigated in a 
forthcoming paper. 

For models having a metallicity equal to or 
greater than 1/10 z (Models 0, 1, 2 and 3 in Ta- 
ble [2), we can separate species into three different 
categories: (a) species with relative abundances 
varying linearly with the metallicity changes e.g. 
CO, H 2 0, CS, SO (see Fig. [TJ and C0 2 , S0 2) 
OCS and H 2 S (see Table [4]); (b) molecules such as 
CN, OH, HCN and HNC showing chemical abun- 
dances that are rather insensitive to the metal- 
licity changes; and (c) species whose abundances 
changes are inversely dependent on z, such as C 2 , 
C 2 H (see Fig. [TJ, and H 2 CO, CH 2 CO and H 2 CS 
(see Table [§. 

Such varieties of behavior have been found in 
the earliest studies of chemical sensitivity to pa- 
rameters ( Pickles fc Williams||1981 1 and in recent 
work ( Bayet et al.|2008b I and it is reasonably well 
understood. Category (a) species may be lost to 
processes independent of metallicity (e.g. CO is 
lost to dissociative ionization with He + ) 

while category (b) species are both formed and 
destroyed in exchange reactions whose reactants 
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Table 1: Standard model parameters (Bell et al. 2006 20071 and references therein. These parameters are 



those used in Models 0-12 except if a different value is specified in Table [2] o r described in the text. 

Grain size 0.1 /xm 

Grain albedo 0.7 

Microturbulence velocity 1.5 kms -1 

Mean photon scattering by grain 0.9 

External FUV radiation intensity (=1©) 1 Habing" 

Cosmic ray ionization rate (=C©) 5.0xl0 _17 s _1 

H2 formation rate coefficient (=R Q ) 3x 10 -18 vTexp( j^q) cm 3 s~ 1 

Gas:dust mass ratio (=d Q ) 100 

Metallicity (=z Q ) solar values b 

Si/H 8.21 xl0~ 7 

Fe/H 3.60 xlO" 7 

Cl/H 1.10 xlO" 7 

Na/H 8.84 xl0~ 7 

Ca/H 5.72 xlO" 10 

a : The unit o f the standard Interstellar Radiation Field (ISRF) intensity is I© = 1.6x10" 
Habingj ( 1968 1; : z = 1 = z corresponds to solar values of the initial elemental abundance ratios (see 



3 erg cm 2 s 



Table [3]) while z = 1/10 z Q means that the solar values of the initial elemental abundance ratios have been 
all divided by the same factor (of 10 in this example). 



depend on metallicity and are therefore largely 
neutral to metallicity changes (For example, OH 
depends mainly on oxygen atoms for its formation 
and on exchange reactions with O and other atoms 
for its loss). Formation schemes for some category 
(c) species depend on the production of C + ions 
from CO, a molecule proportional to z, but loss 
by a succession of exchange reactions resulting in 
a strong inverse dependence on z; thus, there is 
an overall inverse dependence of the abundance of 
such species on z. 

4.2. Initial elemental abundance ratios 

We have investigated the consequences of mak- 
ing changes in the initial elemental abundance ra- 
tios using Models 9, 10 and 11 (see Table [2]). These 
three models correspond to a study of the influ- 
ence of a decrease of nitrogen (Model 9), and an 
increase of both the oxygen and sulfur content 
(Model 10 and 11). The values of the correspond- 
ing initial elemental abundance ratios are listed in 
Table [31 Model 9 refers to values from ChicfR & 



Limongi (20021 while Model 10 and 11 used the 



values suggested by Heger & Woosley ( 2002[ ) and 



Umeda & Nomoto (20021, respectively. Several 



time-dependent chemical trends are shown in Fig. 

m 



Neither oxygen- and sulfur-rich environments 
nor nitrogen-poor environments strongly affect the 
values of the abundances of CO, OH, H 2 S, HCO + , 
H2O and H30 + . Indeed, after 10 5 yr, the scatter in 
their relative abundances is small, within a factor 
of < 10. The same behavior is seen for the CO2 
relative abundance. These species are thus likely 
to be irrelevant as tracers of initial metal gas con- 
tent. On the contrary, SO and CS (see Fig. J2I as 
well as OCS and SO2 abundances (see TablefSTare 
maximum for the oxygen- and sulfur-rich models 
and the variation of their relative abundances is 
important (within a factor of < 100). This makes 
them very sensitive tracers. For Model 0, it is 
shown that CN, HNC, HCN, C 2 and C 2 H reach 
their highest relative abundances having in addi- 
tion significant variations (factor > 10 3 ) as com- 
pared to Model 9 (nitrogen-poor model). These 
last two species are seen particularly sensitive to 
the nitrogen gas content since they are likely to 
only survive in nitrogen-rich environment (for a 
N/H ratio higher than 6.5 xl0~ 5 ). The models 
having the highest O/H and S/H ratios do not 
reveal H2CO, H2CS and CH2CO as particularly 
good tracers. Their corresponding relative abun- 
dances reach their lowest values for Models 10 and 
11 while their maximum correspond to Model 0. 
The decrease in relative abundance seen for some 
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10* 10 s 10 s 
time {yrs) 

Fig. 1. — Influence of metallicity variations on 
some molecular abundances. We plot in a loga- 
rithmic scale the relative abundances n(X)/n# of 
various species, X, with respect to the time (in 
yrs) . Results from the first five models in Table [2] 
(Models to 4) are plotted in this figure: Model 
(solar metallicity) is represented by black lines, 
Model 1 (1/2 z ) by dash- two-dots lines while 
Models 2 (1/4 z ), 3 (1/10 z Q ) and 4 (1/100 z Q ) 
are symbolized by dash-one-dot lines, dashed lines 
and dotted lines, respectively, showing the effect 
of decreasing the metallicity. 



1 D 1Q 



of them when the O/H and S/H elemental abun- 
dance ratios are increased, is so pronounced (e.g. 
up to a factor of 1000 for H 2 CO) that the abun- 
dances fall under the limit of detectability. These 
molecules are thus likely to be undetectable in 
oxygen- and sulfur-rich primordial gas while po- 
tentially observable in nitrogen-rich environments. 
The trends we have identified are summarized in 
Table 

4.3. A„ 

As expected, and shown on Fig. [3j the abun- 
dances of many molecular species are much smaller 
at an A„ = 3 (Model 5 in Table |2| than at A„ = 8 
(Model in Table[2|. For example, H 2 0, HNC and 
SO molecules show fractional abundances smaller 
by more than two orders of magnitude at A„ = 
3 than at A„ = 8. For CS and S0 2 , the reduc- 
tion factors are even larger (see Fig. [3] and Ta- 
ble |6j. These species are very sensitive to A„ 
and good tracers of opacity. At A v =3, how- 
ever, some species become undetectable since their 
corresponding relative abundance drop below the 
probable limit of detectability. 

The strong response to A„ is due to the ability 
of FUV photons to penetrate the translucent cloud 
conditions (A v w 3) but not much deeper. CO is 
self-shielding and generally insensitive for A v > 
1. Some other molecules (C 2 , C 2 H) depend on 
dissociative ionization of CO, and so are also fairly 
insensitive to A v . 

4.4. Cosmic ray ionisation rate 

Using Table [7] and Fig. [4] we identify five 
groups of species showing various degrees of sen- 
sitivity to the change in the cosmic ray ionisation 
rate (from Model with C = 5 x lO^^s" 1 to 
Model 6 with C = 5 x 10 -15 s-\ see Table [2). The 
first group, including OCS, S0 2 or H 2 CS, contains 
the most sensitive species. Indeed these molecules 
have their respective chemical abundances reduced 
by a factor of ~ 10 4 when £ increases by two or- 
ders of magnitude. This approximately quadratic 
response to the enhancement in £ make them very 
good tracers of cosmic rays, especially in very ac- 
tive star-forming regions or where an AGN contri- 
bution occurs. However, this very high sensitivity 
leads certain species to become undetectable at 
high £. This is the case for OCS and S0 2 whose 
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Fig. 3. — Influence of variations on some molecular abundances. We plot in a logarithmic scale the 
relative abundances n(X)/n# of various species, X, with respect to the time (in yrs). The upper two panels 
(dashed lines) are dedicated to Model 5 (A„ — 3 in Table [2]) while the lower two plots (black lines) present 
the results of the solar model (Model in Table [2]) for A v = 8. 



relative abundances thus fall below our (arbitrary) 
limit of detectability of 1 x 10~ 12 for models with 
C= 5 xlO-^s" 1 . 

A second category contains species such as SO, 
H 2 S, CH 2 CO, CS, HNC or HCN whose abun- 
dances respond approximately linearly to changes 
in £. These species are good tracers of cosmic rays 
as long as their fractional abundances are above 
the detectability limit. 

A third category of species includes molecules 
such as CO or H2O which show modest reduction 
of their chemical abundances to increases in (. 

Species such as CN, OH or HCO + are fairly 
insensitive to changes in £. 

Finally, we identify species whose fractional 
abundance increases slightly (by a factor of ~ 10; 
such as H3O" 1 ") or more heavily (by a factor ~ 100; 
such as C2 and C2H) in response to an increase of 
C- Ratios of abundances of molecules from the first 
and fifth groups would be highly sensitive tracers 
ofC 

The higher value of £ used here generates loss 
rates for most species that exceed other loss rates, 
and hence most abundances decrease for higher (. 
However, the larger £ also tends to increase the 



rate of forming ions such as C + , so that transient 
species such as C2 and C2H increase. 

4.5. FUV radiation field 

As expected, increasing the FUV radiation field 
by two orders of magnitude leads to a general 
decrease of molecular abundances (see Fig. [5] 
comparing Models and 7 listed in Table [2] see 
also Table [8]) . Sulfur-bearing species appear also 
very sensitive to the changes in the radiation 
field which makes them very good tracers of FUV 
flux. Molecules such as H2S or H2CS show rela- 
tive abundances linearly reduced with respect to 
the FUV increase while the CS, OCS and S0 2 
molecules are more heavily reduced, by a factor of 
~ 1000. However, SO is less sensitive. Molecules 
CO, CN, OH, HNC, HCN, HCO+ and H 3 0+ seem 
to be rather insensitive to the changes in the FUV 
radiation field. 

The small increase of the relative abundances 
of molecules such as C2 and C2H is due to the in- 
jection of free carbon from the photodestruction of 
carbon-bearing species. This injection stimulates 
hydrocarbon growth. 



G 




10 ! 10* 10 s 1CJ" 10'° 10 1 10* IC S 10* 10 
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Fig. 4. — Influence of cosmic ray ionisation rate variations on some molecular abundances. We plot in a 
logarithmic scale the relative abundances n(X)/n# of various species, X, with respect to the time (in yrs). 
The upper two panels (dashed lines) are dedicated to Model 6 (C enhanced 100 times) while the lower two 
plots (black lines) show results for Model 0. 




10 2 10* 10 1 10" 10 ,B 10 z 10* 10 e 10* 10 
time (yra) time {yfs} 



Fig. 5. — Influence of FUV radiation field variations on some molecular abundances. We plot in a logarithmic 
scale the relative abundances n(X)/n# of various species, X, with respect to the time (in yrs). The upper 
two panels (dashed lines) are dedicated to Model 7 (FUV radiation field enhanced 100 times) while the lower 
two plots (black lines) show results for Model 0. 
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Fig. 2. — Influence of initial elemental abundance 
ratios variations on some molecular abundances. 
We plot in a logarithmic scale the relative abun- 
dances n(X)/n# of various species, X, with respect 
to the time (in yrs). Only the solar (Model 0) and 
the last three models from Table [2] (Models 9 to 
11) are plotted. Model is represented by black 
lines, Model 9 by dotted lines while Models 10 and 
11 are symbolized by dash-one-dot lines, dash-two- 
dots lines, respectively. The values of the initial el- 
emental abundance ratios for each model are given 
in Table H 



4.6. Density 

We can roughly separate species into four main 
categories, depending on their behavior when the 
density changes. Using Fig|6] and Table [9j we 
identify a first category including species rather 
insensitive to gas density changes. This includes 
molecules such as CO, H 2 0, H 2 CS, CS, HCO+, 
H 3 + , HCN or HNC whose abundance variations 
are below a factor of ~5. This category is fol- 
lowed by the molecules SO and OH which show 
moderate change in their chemical abundances by 
a factor of « 10 (especially after 10 6 yr, see Fig. [fjj. 
A third category includes species that respond to 
an increase in gas number density from 10 4 to 10 6 
cm -3 with heavy reduction (CN and C2, factors of 
100), and very heavy reduction of their fractional 
abundances (C2H, CH 2 CO or H 2 CO, factors from 
>100). These very heavy reductions usually lead 
to species becoming non-detectable (see Table |9| . 
For this reason, C 2 H, CH 2 CO H 2 CO cannot be 
considered as good tracers of density. On the 
contrary, CN and C 2 molecules could be seen as 
very good tracers of gas density since they are lin- 
early linked with density changes. A last group of 
species able to be recognized as very good trac- 
ers of density, contains the molecules C0 2 , OCS, 
H 2 S and S0 2 . The relative abundances of these 
species indeed vary linearly (factor of 100) with 
the change in gas density, except for the S0 2 which 
varies quadratically with the gas density, promot- 
ing this molecule to the best (detectable) tracer of 
n(H 2 ). 

However these results take no account of the 
fact that freeze-out of gas phase species has been 
omitted from these models. For gas of num- 
ber density n# and metallicity z, the freeze-out 
timescale is sa 10 6 yr(10 4 /n//)(z Q /z). Thus, at 
n# = 10 cm -3 and z—zq the timescale is around 
one million years, a plausible lifetime for gas in a 
PDR. However, at n# = 10 6 cm~ 3 , the timescale 
becomes unfeasibly short unless the metallicity is 
very low. Thus, unless desorption processes are 



operating efficiently (Roberts et al. 20071 signif- 



icant freeze-out will be occurring in gas at very 
high densities. 
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Fig. 6. — Influence of density variations on some molecular abundances. We plot in a logarithmic scale the 
relative abundances n(X)/n# of various species, X, with respect to the time (in yrs). The upper two panels 
(dashed lines) are dedicated to Model 8 (number density of 10 6 cm" 3 ) while the lower two plots (black lines) 
show results for Model 0. 



5. Predicted PDR tracers for selected 
galaxy types 

In this section, we have selected several fairly 
well-known galaxies as illustrating examples of the 
starburst, SB+AGN and high-z galaxies. These 
galaxies are compared with our selection of mod- 
els in an attempt to identify, for each case, the 
best molecular tracers of dense PDR gas. This 
study is of particular interest for the preparations 
of the future observational programs of Herschel 
and ALMA. 

An interesting type of galaxy concerns objects 
such as M 83, NGC 253, i.e. typical starbursts. 
It is commonly accepted that one may reproduce 
such environments by using a FUV radiation field 
enhanced as compared to normal spiral galaxies 



( |Israel k Baas||2001| |Bayet et alT]|2006[ |2008b[ ) 



M 83 and NGC 253 show solar metallicity (see 
Origlia et al.|2004| |Smith et al.|2006"l and |Zaritsky| 



et al. 



1994 



respectively), leading to a choice of 
standard elemental initial abundance ratios and 
density. Model 7 appears consequently a good 
model likely to represent this morphological type 
of galaxy. 

Another morphological type of galaxy that we 



have investigated concerns galaxies hosting an ad- 
ditional source such as active nucleus (SB+AGN). 
Objects such as the galaxies Arp 220 ( |Soifer et"aL 



1984 



1986 



Joy et al. 1986) or NGC 3079 (Ford et al. 



Sosa-Brito et al. 2001 ) are well-known exam- 



ples of this category. Indeed, whether Arp 220 has 
an AGN or not is still a very highly debated topic 



(see for instance 


Downes k Eckart|2007 Sakamoto 


et al. 2008 


Aalto et al. 2008 


1. Although NGC 3079 



certainly contains an AGN, it seems that also 
for this source the more extended SB may domi- 
nate the molecular gas properties (see for instance 
Kohno et al. 2007). Following the parameter 



choices from Bayet et al. ( 2008b[ ), we have adopted 
here Model 6 to represent these environments. In- 
deed, in this model, the cosmic ray ionisation rate 
is enhanced 100 times (( = 5 x 10~ 15 s _1 ) as com- 
pared to the standard values while the FUV ra- 
diation has been maintained at 10 3 I©. The se- 
lected metallicity is solar. The density, the ini- 
tial elemental abundance ratios and the A v values 
are standard. The model adopted here is thus an 
attempt to simulate qualitatively XDR-likc envi- 
ronments. To reproduce "pure" AGN sources, full 



XDR models (|Meijerink et al.||2005| |Meijerink fc 
Spaans 2005 Spaans k Meijerink 2005) are how- 
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ever required. 

Finally, we want to investigate galaxies at high 
redshift. Typical examples to be considered are 
the Cloverleaf QSO and APM 08279 with then- 
respective redshifts of z r ed. ~ 2.6 and z r ed. ~ 
3.9. For these sources, the physical properties 
are unknown or highly uncertain, even if sev- 



eral molecules have already been detected (Garcia- 



2007 Wcifi ct al. 2007 Riechers et al. 2009). The 



Burillo et al. |[2006l [Guelin et al.||2007[|Krips et al. 



presence of both dense and less dense gas make 
these sources particularly good illustrating ex- 
amples for our study. Some of these measure- 



ments ( 


Gao & Solomon 


2004 


Wu et al. 


2005 


Wagg et al.|2005 


1 suggest that star formation may 



be active, even if these galaxies are affected by 
an AGN. |Garcia-Burillo et al. (2006) have sug- 



gested that the unusual chemistry observed may 
be due to a combination of star formation and 
AGN (PDR/XDR) chemistries. In an attempt to 
describe such sources, we have run an additional 
model presented in Table [2] (Model 12). This 
model shows a high value of the cosmic ray ion- 
ization rate, a high value of the radiation field in- 
tensity, and a relatively low metallicity (all with 
respect to the Milky Way). 

We summarize in Table |10| the species predicted 
to be the best molecular tracers of dense PDR gas 
in each case (above the limit of detectability arbi- 
trarily taken to be equivalent to a fractional abun- 
dance of lxlO -12 ). Here our goal was to give the 
reader some observational clues derived from this 
modelling work. Our choice of galaxies are simply 
examples; none of the models we are presenting in 
Table [10] are meant as detailed studies of own in- 
dividual galaxy. Together with the overall trends 
listed in Sect. [4j Table [10| should enable a more 
detailed description of the physical properties of 
the galaxy to be made. 

Whichever model (6, 7 or 12) is used, CO and 
H 2 molecules are predicted to have high frac- 
tional abundances (from « 3 x 10~ 8 up to 10 -4 ). 
More surprisingly, for Models 6 and 7, CS, CN 
and OH also show relatively high fractional abun- 
dances, ranging between 10 -9 and 10 -8 . Ac- 



cording to the estimates of Lintott et al. (20051 
these abundances are sufficiently large that unre- 
solved active galaxies should be detectable in these 
species. Observational tests of the predictions 
listed in Table [10] have been already performed by 



Bayet et al. (2008a I. They have revealed the im- 
portance of CS fines in starburst galaxies. Further 
observational tests of predictions made in Table [T0| 
will be detailed in forthcoming papers. For Model 
12, CN and CS have chemical abundances which 
drop severely by at least two orders of magnitude 
while the OH relative abundance increases by a 
factor of 10 relative to Model 0. This feature again 
reinforces the interest of detecting these species, 
especially at high redshift. 

Finally, it is interesting to note that molecules 
such as HCN, HNC and HCO+ do not show es- 
pecially high relative abundances in the models 
investigated in this section. Indeed, the HCN and 
HNC chemical abundances for Models 6, 7 and 
12 fall to « 10" u . Moreover, the HCO+ abun- 
dance varies from 2.0 xl0~ n (Model 7) to 4.0 
xl0~ 10 (Model 12). However, we note that us- 



ing a full XDR model ( 


Meijerink & Spaans 


2005 


Spaans & Meijerink |2005 


Meijerink et al. 


20071 



might change the results of Model 6. Indeed, in 
a full XDR model the cosmic rays considered are 
much more energetic than the ones normally used 
in PDR models. 

6. Discussion and conclusions 

The main results of this study are the trends 
in the chemistry in response to changes in the 
physical conditions (model parameters) summa- 
rized in Tables [4] to [9j In addition to the evo- 
lution of the chemistry, Table [10] may be impor- 
tant in the context of the preparation of future ob- 
servational programs (Herschel, ALMA). Our re- 
sults focus mainly on the chemistry taking place 
at A^ 8, i.e. in the densest part of the PDR. 
Analysis and discussion on the chemistry taking 
place in the outer region of the cloud (A^ <C 3) 
where the gas temperature and the FUV radiation 
field are higher is beyond the scope of this paper. 



While previous studies (Israel et al. 


1995 


Gerin 


& Phillips 


2000 |Israel & Baas |2001 


20021 2003| 


Bayet et al.|2004 


2006) focused their 


attention on 



determining the physical conditions of a few extra- 
galactic star-forming or active regions using vari- 
ous PDR codes and the emissions of only a couple 
of molecules (mainly neutral carbon and carbon 
monoxide), we provide here a wide and complete 
analysis of the chemistry, and report in details on 
19 species in a large range of parameter space. 
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Our model predictions must be compared with 



observations (e.g. 


Radford et al. 1991 Solomon 


et al.|1992||2004| 


Martin et al. | 


2005 Evans et al. 


2006| Gracia-Carpio et al.|2006| Martin et al.|2006 


Aalto et al. 2007 Imanishi et al. 2007| Knud- 


sen et al. 2007 1 


3 apadopoulos 


20071 to be val- 



idated. Recent extragalactic detections of HCN 
and HCO + in 12 nearby galaxies summarized in 



Krips et al. (2008) provide useful estimations of 
molecular abundances. When the source hosts an 
AGN contribution, the relative HCN abundance 
they predicted from an LVG analysis of their data, 
is in agreement with our predictions (model 6), to 
within a factor of 10. For starburst-dominated 
galaxies, we also obtain a good agreement. As an- 
other exampl e, H 3 Q + has been rece ntly detected 
in Arp 220 ( |van der Tak et al.||2008) . The H 3 0+ 
abundance in this SB+AGN source is estimated 
to be 2-lOxlO -9 . Model 6 predicts an abund ance 
of H30 + lower by a factor of 10 



However, 



der Tak et al. ( 2008 1 mentioned that their abun- 
dance estimation is somewhat hampered by source 
size uncertainties, which we think may explain the 
discrepancy between our predictions and their ob- 
servations. We conclude that our models appear 
to be in harmony with the limited observational 
data. 

It is also particularly useful to compare the 
molecular tracers predicted by dense star-forming 



core models (IBayet et al. 2008b I with those pre 



sented in this article (derived from the UCL_PDR 
code) . Indeed, it is essential for better understand- 
ing and diagnosing the nuclear energy source and 
activity in galaxies to be able to disentangle the 
PDR emission from that produced by dense star- 
forming cores. This should lead to a better de- 
termination of the excitation conditions in such 
regions. These studies are of importance because 



they concern not only the local Universe (Kohno 



2008 



et al : ||2001[|Gracia-Carpio et al.||2006| |Baan et aT 



Krips et al. 20081 but also could help us 



to better understand the galaxy formation in the 



Universe at higher redshifts ( 


Ivison et al.|2005b|a 


2008 


Younger et al. 


2008b|a 


). Comparing model 



predictions we thus have obtained that for the 
starburst galaxy case, S0 2 and H 2 S are predicted 
to be likely inappropriate tracers of PDR regions 
while they are predicted as good tracers of dense 
star- forming cores. By contrast, HCO + is pre- 
dicted to be undetectable in dense star-forming 



cores while enhanced in PDR-dominated galaxies. 
For high redshift sources, H2CS and H2CO are 
both predicted to be inadequate molecular tracers 
of PDRs while both are good signatures of dense 
star- forming cores. 

By comparing the model predictions for star- 



burst galaxies (Bayet et al. 2008b I and for the 



present study of PDR chemistry, we note that 
while the molecular component of AGN-excited 
galaxies is typical of PDRs, galaxies dominated 
by clusters of massive stars (SB) contain molec- 
ular tracers of both PDR and dense star-forming 
gas. Thus, in principle, it should be possible to 
distinguish between SB+AGN and starburst exci- 
tations. 

In summary, we have explored PDR chemistry 
for a range of parameters representing galaxies 
with intense FUV or cosmic ray sources. Our main 
results demonstrate the sensitivity of the chem- 
istry to the local physical conditions, and show 
that molecular observations can be used to de- 
termine those conditions. We have adopted very 
crude estimations of the relevant physical condi- 
tions in several galaxy types and use our models 
to predict molecular species that may be the best 
tracers of those conditions. 
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Table 2: Input parameters of the UCL_PDR models used to perform this study (see Sects. [3] and [3]). The 
abbreviation "ST" represents the standard values while the abbreviations "CL02", "HW02" and "UN02" 
are the initial elemental abundance ratios references detailed in Table [3] This table does not present all the 
input parameters of the UCLJPDR code for each model but only lists the parameters set to values different 
from the standard ones. 



Model 


Metallicity 


Gas-to-dust 


H2 form. 


A„ Ini. Elem. 


FUV rad. field 


c 


Gas Density 




(-0) 


mass ratio (d©) 


rate coeff.(R Q ) 


Abund. ratios 


(Io) 


(Co) 


(cm -3 ) 





1 


1 


1 


8 ST 


10 3 


1 


10 4 


1 


1/2 


1/2 


1/2 


8 ST/2 


10 3 


1 


10 4 


2 


1/4 


1/4 


1/4 


8 ST/4 


10 3 


1 


10 4 


3 


1/10 


1/10 


1/10 


8 ST/10 


10 3 


1 


10 4 


4 


1/100 


1/100 


1/100 


8 ST/100 


10 3 


1 


10 4 


5 


1 


1 


1 


3 ST 


10 3 


1 


10 4 


G 


1 


1 


1 


8 ST 


10 3 


100 


10 4 


7 


1 


1 


1 


8 ST 


10 5 


1 


10 4 


8 


1 


1 


1 


8 ST 


10 3 


1 


10 6 


9 


1 


1 


1 


8 CL02 


10 3 


1 


10 4 


10 


1 


1 


1 


8 HW02 


10 3 


1 


10 4 


11 


1 


1 


1 


8 UN02 


10 3 


1 


10 4 


12 


1/10 


1/10 


1/10 


8 ST/10 


10 5 


100 


10 4 



T able 3: Initial abundance r atios used in Table |2) T he a bbreviations "CL02" "HW 02" and "UN02" refer 



tc 


) Chieffi & Limongi 


(2002 


, Heger & Woosley 


(2002|) and Umeda & Nomoto 


(2002 


1, respectively (see Sect. 


. The standard initial elemental abundance ratios ("ST") are from 


Sembach & Savage ( 


1996 


); Sofia et al. 


<] 


J97 


); 


Meyer et al. 


1998); 


Snow et al. 


(2002); Knauth et al. ( 


2003 


. Models using the "ST" values of the 



initial elemental abundance ratios are the models corresponding to z = z Q . 





ST 




CL02 




HW02 




UN02 




C/H 


1.4x10 


-4 


1.4x10" 


4 


1.4x10" 


4 


1.4x10" 


4 


O/H 


3.2x10 


-4 


4.54x10 


-4 


1.53x10" 


-3 


1.18x10 


-3 


N/H 


6.5x10 


-5 


5.99x10- 


-11 


1.58x10" 


-9 


3.24x10 


-7 


S/H 


1.4x10 


-6 


7.59x10 


-6 


1.54x10" 


-4 


4.06x10 


-5 


He/H 


7.5x10 


-2 


7.5x10" 


2 


7.5x10" 


2 


7.5x10" 


2 


Mg/H 


5.1x10 


-6 


1.83x10 


-5 


1.2x10" 


4 


6.18x10" 


-5 
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Table 4: Trends of molecular fractional abundances, computed for different values of metallicity. When a 
single species is denned as "undetectable" or "not detectable" , it means that its relative abundance is below 
the adopted limit of detectability of 1 x 1CP 12 . 



Molecule Response to metallicity changes 

CO, H2O linear tracers of metallicity 

CS, SO for z>l/lOOz , good linear tracers of metallicity 

CN for z>l/lOOz ; insensitive to metallicity changes 

OH, H3O 4 " most abundant at z=l/lOOz ; generally insensitive to z 

HNC, HCN for z>l/10z Q , insensitive to metallicity changes 

HCO + always above the limit of detectability; insensitive to metallicity changes 

C2, C2H for z>l/lOz , inversely dependent on z changes 

CO2, OCS good linear tracers of metallicity except for z<l/100z Q (undetectable) 

S02,H 2 S for z>l/100z Q , good linear tracer of metallicity 
H2CS, CH2CO for z>l/10z Q , inversely dependent on z changes 

H2CO for z>l/lOOz , inversely dependent on z changes 



Table 5: Trends of molecular fractional abundances, computed for different values of the initial elemental 
abundance ratios. Comparison of the PDR chemistry for various values of initial abundances (see Table [3] 
and Models 9, 10 & 11 in Table [2]) while other parameters have the standard values (see Table [TJ. 



Molecule Response to initial elemental abundance ratios changes 

CO, H2O insensitive to changes of any initial elemental abundance ratios 

CS sensitive to changes in both S/H and O/H; variations within a factor of > 50 

SO very good tracer of S/H; variations up to a factor as 1000; very sensitive to O/H changes 

CN highest abundance for N-rich models; very good tracer of N/H 

OH insensitive to any changes in initial elemental abundance ratios 

HNC, HCN very good tracers of N/H 

HCO + , H30 + poor tracer of initial elemental abundance ratios 

C2, C2H surviving only if N-rich environments 

CO2 insensitive in changes of any initial elemental abundance ratios 

OCS highest abundance for S- and O-rich models; variations within a factor of sa 100; good tracer 

SO2 very good tracer of S/H and O/H 

H2S insensitive to any changes in initial elemental abundance ratios 

H 2 CS bad tracer of S/H and O/H; lowest abundance for the S-and O-rich models 



CH2CO, H2CO abundant when N-enriched environment, poor tracers of S/H and O/H 
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Table 6: Trends of molecular fractional abundances, computed for two values of the opacity (A„). Comparison 
of the PDR chemistry for A„ = 8 (Model 0) and A„ = 3 (Model 5) while other parameters have the standard 
values (see Table [lj . 



Molecule 


Response to opacity changes 


CO 


insensitive to opacity changes 


H 2 0, HNC, HCN, SO 


very good tracers of A„ 


CS 


at A„ = 3, very good tracer of the opacity 


CN, HCO+, H 3 0+ 


not very sensitive to opacity changes 


OH 


insensitive to opacity changes; poor tracer of A„ 


C2, C2H 


insensitive to opacity changes 


C0 2 

OCS, H 2 CS 
S0 2 

H 2 S, CH 2 CO 
H 2 CO 


at A„ = 3, abundance decreases by a factor of 10 4 ; very good tracer of A„ 

at A v = 3, reduced by a factor of 10 4 ; become undetectable 

at A v = 3, reduced by a factor of 10 6 ; become undetectable 

at A v = 3, decreased by a factor >10 2 ; become undetectable 

when A„ = 3, abundance reduction by a factor of 10; become undetectable 



Table 7: Trends of molecular fractional abundances, computed for different values of the cosmic ray ionization 
rate (£). Comparison of the PDR chemistry for C = 1£ = 5 xlO^V 1 (Model 0) and £ = lOOxC© = 5 
xl0~ 15 s _1 (Model 6) while other parameters have the standard values (see Table [lj. 

Molecule Response to cosmic ray ionisation rate changes 

CO, H 2 reduced by a factor of 10 with £; possible tracers of C 

CS, SO, HNC, HCN heavy (linear) abundance reduction by factor of 100 with £; good tracers of £ 

CN, OH, HCO+ insensitive to £ 

H 3 + slight increase by factor of 10 with £ 

C2, C2H heavy (linear) abundance increase by a factor of 100 with £ 

CO2, H2CS very heavy reduction by a factor of ~ 10 4 ; very good tracer of £ 

OCS, SO2 very heavily reduced by factor of ~ 10 4 with £; become undetectable 

H 2 S, CH 2 CO heavy (linear) abundance reduction by factor of 100 with £; become undetectable 

H2CO insensitive to £ 
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Table 8: Trends of molecular fractional abundances, computed for different values of the FUV radiation field. 
Comparison of the PDR chemistry for I = 10 3 I (Model 0) and i' = 100 x 1= 10 5 x I (Model 7) while 
other parameters have the standard values (see Table [ll . 



Molecule Response to FUV radiation field changes 

CO insensitive to FUV changes 

H2O reduced by a factor of 100 (linearly to the FUV change); good tracer 

CS heavy abundance reduction by a factor of 1000 with FUV; very sensitive to FUV changes 

SO moderately reduced by a factor of 10 with FUV; possible tracer of FUV 

CN, OH, HNC, insensitive to the changes in FUV 

HCN, HCO+, H 3 0+ insensitive to the changes in FUV 

C2, C2H slight abundance increase by a factor < 10 with FUV changes 

CO2, H 2 S, H 2 CS abundance decreases (linearly) by a factor of 100 with the FUV changes; good tracers 

OCS, SO2 heavy reduction by a factor of 1000; very sensitive to FUV changes; very good tracer 

CH 2 CO, H 2 CO insensitive to FUV changes 



Table 9: Trends of molecular fractional abundances, computed for different values of the density. Compar- 
ison of the PDR chemistry for n(H 2 )= 10 3 cm -3 (Model 0) and n(H 2 ) = 10 6 cm~ 3 (Model 6) while other 
parameters have the standard values (see Table [I]). 



Molecule Response to density changes 



CO, H2O, CS, HNC, HCN abundance variations within a factor of < 4; poor tracers 

SO abundance increase by a factor of ~10 with density; possible tracer of density 

CN, C 2 decrease by a factor of 100 (linearly inverse with the change of density); very good tracer 

OH decrease by a factor of ~10; possible tracer of density 

HCO + , HaO" 1 " abundance variations within a factor of < 4; insensitive to density changes 

C2H decrease by a factor of > 100; undetectable at high gas density 

CO2, OCS, H 2 S increase by a factor of >100; very good (linear) tracers of n(H 2 ) 

S0 2 increase by a factor of >1000; very good (quadratic) tracer of n(H 2 ) 

H2CS abundance variations within a factor of ~ 2; insensitive to density changes 

CH 2 CO decrease by a factor of >10 5 ; undetectable at high density 

H 2 CO (linear) decrease by a factor of >100; undetectable at n(H 2 )=10 6 cm~ 3 



15 



2009, MNRAS, in preparation 



Table 10: Detectability of 19 species likely to trace 
the dense PDR gas component in three categories 
of galaxies (see Sect. |5|. The limit of detectability 
has been taken to be (n(X)/n#) = 1 xl0~ 12 , as 
is typical for molecules in the Milky Way. Under 
this limit the species are not detectable (symbol 
— ). Otherwise, they are marked with the symbol 
+ . We have separated the species plotted in Figs. 
[1] to [2] from those only listed in Tables [4] to [5] 





Model 7 


Model 6 


Model 12 


Galaxy type 


Starburst 


SB+AGN 


high redshift 


CO 


+ 


+ 


+ 


H 2 


+ 


+ 


+ 


CS 


+ 


+ 


+ 


SO 


+ 


+ 


+ 


CN 


+ 


+ 


+ 


OH 


+ 


+ 


+ 


HNC 


+ 


+ 


+ 


HCN 


+ 


+ 


+ 


HCO+ 


+ 


+ 


+ 


H 3 0+ 


+ 


+ 


+ 


c 2 


+ 


+ 


+ 


C 2 H 


+ 


+ 


+ 


C0 2 


+ 


+ 


+ 


ocs 


+ 






so 2 








H 2 S 








H 2 CS 


+ 


+ 




H 2 CO 


+ 


+ 




CH 2 CO 


+ 
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